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Abstract. This article is the first in the series that analyze about 1.87 million periods of PSR B0531+21 (Crab 
pulsar), observed by the PC A detector aboard the RXTE x-ray observatory. The Crab pulsar's x-ray light curve 
shows little variation over time scales ranging from days to a period (33.46 milli seconds) . The standard deviation 
of its x-ray flux variation is w 0.7% of its mean value, which is negligible compared to its radio flux variations. The 
phase resolved power spectrum of pulse to pulse x-ray flux variation shows no spectral feature; an upper limit to 
the peak of any possible broad spectral feature is 0.06% of the mean power. The x-ray fluxes in the two components 
of its integrated profile are unrelated to each other; their linear correlation coefficient is 0.0004±0.0010. "Giant 
pulses" that are routinely seen at radio wavelengths are absent here. This work sets very strong constraints on 
the connection (if any) between the flux variations at radio and x-ray energies, for example due to variation in the 
degree of coherence of the basic emitters. Its phase resolved x-ray flux variation shows a weak correlation with the 
integrated profile. If confirmed, this might be an important clue to understanding the x-ray emission mechanism 
of Crab pulsar. 
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1. Introduction 

The intrinsic temporal flux variations of pulsars are prob- 
ably an important clue to understanding their emission 
mechanism at any wavelength. In fact, the flux variations 
at high energies (optical, x-ray, 7-ray) might provide im- 
portant constraints to the emission mechanism at radio 
wavelengths also (Cheng et al, 1986b| Kawai et al 1991 



Lundgren et al |1995|; M offet & Hankins |1996[ Patt el al 



1999). Patt el al (1999) searched for period to period flux 



variations in about 105 000 periods of Crab pulsar data 
at x-ray energies (in the range 1 to 10 kilo electron volts 
(KeV)), with a time resolution of 100 micro seconds {/is), 
obtained by the PCA detector aboard the RXTE x-ray 
observatory. This work reports the results of analyzing 
1868112 periods of Crab pulsar from the same instru- 
ment, with a time resolution of 3.815 /zs, in the energy 
range 13.3 to 58.4 KeV. 

The RXTE data archive was searched using the XDF 
tool, for public data acquired by the Proportional Counter 
Array (PCA). A uniform set of 23 data files were found ob- 
served during August/September 1996, their Obsld num- 
bers ranging from 10203-01-01-00 to 10203-01-03-01. They 
were obtained in the EVENT mode (XTE_SE), combining 
events from all five Proportional Counter Units (D[0~4]), 
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and also from both halves of all three Xenon anode layers 
of eachPCU (X1I/X1R"X2L"X2R"X3I/X3R). Channels 
50 to 249 of the PCA were also combined, which corre- 
sponds to the energy range 13.3 to 58.4 KeV. 

The first phase of data analysis used the FTOOLS soft- 
ware. First, the Good Time Intervals (GTI) were obtained 
for each data file by using the MAKETIME tool on the 
corresponding XTE filter file; the selection criterion were 
(a) pointing OFFSET less than 0.02°, (b) elevation (ELV) 
greater than 10°, and (c) all five PCUs to be switched on. 
Next, the GTI extension of each data file was edited to in- 
sert the above GTI values. Then the FSELECT tool was 
run to filter out data outside these time ranges. Next the 
SEFILTER tool was run with the M[1]{1} option (with- 
out bypassing the FSELECT tool) to retain only the valid 
photon events. Then the FXBARY tool was run using the 
orbit file for that day, to convert the arrival times of pho- 
tons from the Terrestrial Time system (TT) to the so- 
lar system barycenter system (TDB). Penultimately, the 
SEEXTRCT tool was used to obtain the light curve for 
each file, in time intervals of 1.010895 milli second (ms), 
which is 265 times the basic time resolution of the data. 
Finally, the PCADTLC tool was used to correct the light 
curve for dead time of the PCA; before this the corre- 
sponding Standard 1 files were also converted to the TDB 
system for consistency. 
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Fig. 1. Integrated profile of Crab pulsar after summing 
1868112 periods from 23 files. The abscissa is time (also 
called phase) within the period (in ms), while the ordinate 
is the average number of photons obtained in one synthe- 
sized time sample (1.013967 ms). The mis-alignment of 
the profiles from file to file does not exceed half a synthe- 
sized time sample. The dotted curves represent the two 
peaks modeled as Gaussian. 



Time samples having incomplete exposure were 
deleted. These occurred naturally at the beginning and 
end of each light curve, and also whenever the RXTE ob- 
servatory shut off some PCUs, for technical reasons. Each 
light curve was then converted into the ASCII format for 
further processing. 

The second phase of data analysis used self-developed 
software. First the power spectrum of each light curve was 
computed to obtain the period of Crab pulsar in that file. 
This was used to separate the light curve into individ- 
ual periods (also called single pulses). Each period has 33 
time samples (also called bins), giving a synthesized sam- 
pling interval that is different from file to file, but is sa 
1.013967 ms. The above separation was done such that 
the photon counts in each original sampling interval were 
not split across more than one synthesized sampling inter- 
vals; otherwise the Poisson statistics of the data would be 
distorted, and would cause problems for some studies as 
discussed later on. The separation into individual periods 
for radio data is much simpler, since standard resampling 
techniques can be used. 

The appropriate period for each data file was ascer- 
tained more accurately by checking for "drift" of the in- 
tegrated profile between the first and second halves of the 
light curve; this part was done iteratively (see Vivekanand 
et al 1998 for details). A straight line fit, to the starting 
epoch (TDB system) of each data file versus the Crab pul- 
sar period in that file, gives a period derivative of 4.208 
(± 0.005) xl0~ 13 s/s, which compares excellently with 
the actual value for Crab pulsar. The standard deviation 
of the periods about the fitted straight line is w 2 nano 
seconds (ns), which is consistent with the expected value; 
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Fig. 2. X-ray light curve of Crab pulsar over 26 days dur- 
ing August/September 1996. The first two plots show the 
average number of photons in the off-pulse window (OFF) 
and on-pulse window (ON) defined in Fig. [|. The last plot 
is the ratio ON/OFF. 



the data files contain typically 50 000 to 98 000 periods, 
and one can recognize a relative shift of a small fraction 
of the time sample between the integrated profiles of the 
first and second halves of a data file. However the periods 
differ systematically by ss 8 ns from those obtained using 
Crab pulsar's ephemeris. It is not clear to this author why 
this systematic difference should occur, but this does not 
affect the rest of the analysis. 

Fig. |] shows the integrated profile of Crab pulsar for 
1868112 periods. The integrated profiles of all files are 
aligned, correct to half a time sample, so that one can an- 
alyze them as if all 23 files have been obtained "in phase" . 
However it does smear the integrated profile to a maxi- 
mum of half a synthesized time sample. Samples 5 to 30 
are considered to represent the on-pulse window, and the 
rest of the seven samples the off-pulse window, although 
the Crab pulsar might emit x-rays all through its period. 

Details of the analysis in the coming sections can be 



found in Vivekanand & Joshi 1997, Vivekanand et al 1998 



and Vivekanand 200C; they will be described only briefly 
in this article. 



2. X-ray light curve 

Figure |^ shows the average number of photons in the off- 
and on- pulse windows of Fig. [fl Each point in the figure 
represents the average number of photons in 5 000 periods 
(ss 167 s); the total duration in the figure is 26.1 days. 
Both ON and OFF counts show correlated hourly as well 
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as daily variation. Assuming that this is most probably 
an instrumental effect, the ratio ON/OFF in Fig. |] would 
represent the true variation of the Crab pulsar x-ray flux. 
The mean value of the ratio is 4.432 while its standard 
deviation is 0.035; so the x-ray flux of Crab pulsar fluctu- 
ates with a standard deviation of 0.035/4.432 w 0.8% of 
its mean value, which is a negligible quantity. This sup- 
ports the reported Crab pulsar behavior at 0.015 to 0.130 



mega electron volt (MeV) energies (Ubertini et al 1994), 
but is contrary to its behavior at w 200 MeV reported 
by Nolan et al ( 1993 ) using CGRO observations, where 
they found the 7 ray flux of Crab pulsar to vary over time 
scales of a month. However the above result (0.8%) would 
be an under estimate if there is significant emission from 
the Crab pulsar during the off-pulse window. 

It would have been interesting to study the existence, 
at x-rays, of the reported very long time scale (of ~ 13.5 
years) variation of the ratio of t he flux es in the two peaks 
(of Fig. 0) at 7-rays (Nolan et al 1993 ); unfortunately this 
can not be done due to the limited data span available 
here. 

3. Pulse to pulse flux variations 

The last section discussed the flux variations of Crab pul- 
sar over time scales of hours and days. The current section 
discusses flux variations from period to period. 

3.1. Spectrum of flux variations 

Figure |^ shows the so called x-ray fluctuation spectrum 
of Crab pulsar. At each of the 33 samples (phases) of the 
integrated profile of Fig. [j], a time series was formed com- 
prising of the x-ray flux as a function of the period number 
in the data file. This was Fourier transformed in arrays of 
length 128 x 1024 periods. The data were centered in 
the array and zero padded, and then a Hamming window 
was applied. To remove long term variations, the data in 
blocks of 32 768 periods were normalized with the mean 



value of this block (see Ritchings 1976 and Vivekanand 
& Joshi 1997 for details). Fourteen data files were chosen 
that had at least 75 000 periods each, totalling to 1 347 028 
periods. Figure ^| shows the power spectrum averaged over 
the 33 spectra, after normalizing each spectrum with the 
variance of its time series. A polynomial of the form 



y = ao + a\x + aix 



a 3 x 



(1) 



was fit to the power spectrum in Fig. ||. the coefficients are 
a = 0.999±0.001, a x = 0.003±0.013, a 2 = -0.003±0.059 
and (13 = —0.004 ± 0.077. The standard deviation of the 
power spectrum with respect to the above fitted curve is 
0.046, which is mainly determined by photon noise. It is 
clear that none of the coefficients are significant except 
the first. 

The upper limit to any possible weak and broad spec- 
tral feature that might be hidden in the data can be com- 
puted to be w 0.06% of the total power in the spectrum. 
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Fig. 3. Average power spectrum of pulse to pulse x-ray 
flux variations of Crab pulsar. The time series, consisting 
of x-ray photons as a function of the period number, was 
Fourier transformed at each of the 33 phases of the inte- 
grated profile; the modulus squared of the 33 transforms 
was averaged. The abscissa in is units of inverse periods, 
extending up to 0.5 due to the Nyquist criterion. A small 
range of frequencies has been removed around the sam- 
pling spectral feature (w 1/10 periods), and two of its 
harmonics. 

Therefore Fig. || is consistent with the Crab pulsar having 
no spectral feature in its pulse to pulse x-ray flux varia- 
tions. 



3.2. Modulation index 

The radio flux of rotation powered pulsars varies signifi- 
cantly from pulse to pulse; this could be due to intrinsic 
flux variation of the sub pulse, as well as random position 
of the sub pulse within the on-pulse window (here one is 
ignoring the flux variations due to propagation in the in- 
terstellar medium). This is characterized by the so called 
modulation index fi, defined as 



M = 



1L 

'J) 



(2) 



where (/) and 07 are the mean and the standard devia- 
tion of the pulsar flux (see Manchester & Taylor 1977); \i 
represents the fractional flux variation of the pulsar, and 
is usually greater than 1.0 at radio wavelengths. 

Figure ^ shows the square of which is the natural 
quantity to average, for all 1868112 periods. First, the 
mean (J) and the variance o 1 are computed at each sample 
in the integrated profile ((/) is plotted in Fig. [I] and the 
top panel of Fig. |i|). Now a 2 has contribution from two 
sources: 



K 



= (!)+< 



(3) 



where a\ is the variance due to Poisson statistics of pho- 
tons, and equals the mean number of photons (I), while er| 
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Fig. 4. Square of the estimated modulation index (/z 2 ) of 
the Crab pulsar x-ray flux, at each sample of the inte- 
grated profile. The top panel is the same as Fig. [l| The 
vertical bars represent two standard deviation error bars. 
Since depends upon the small difference between two 
much larger quantities (eq. 3), it can be negative also 
depending upon measurement errors. The distribution of 
these negative values about the value 0.0 is consistent with 
their rms errors, as expected. 



represents the fluctuation of the average intensity of the 
pulsar; the two are referred to as photon noise and wave 
noise, respectively (Goodman f 985 ). The modulation in- 
dex squared in Fig. |J was estimated by subtracting the 
mean intensity (I) from the estimated variance a 2 , and 
then dividing by (I) 2 at each sample of the integrated 
profile. Because of the PCA dead time correction to the 
data, the photon rate of Crab pulsar in any time sample 
is about 6% higher than the corresponding integer value 
of photons. This correction was estimated self-consistently 
by averaging a 2 — (I) over all 33 samples of the integrated 
profile, for each file. The average value of the correction for 
all 23 files is 1.062 with standard deviation of 0.003. The 
mean flux (J) at each sample was scaled by this constant 
(for that data file) before subtracting from a 2 . This con- 
stant is not dependent upon the average flux (I) at each 
sample, so the above procedure is unlikely to introduce 
artifacts in the /i 2 of Fig. [|. 

The average value of [i 2 in Fig. | is -0.0002 ± 0.0040, 
which is is consistent with zero as expected. The x 2 °f 
/i 2 of Fig. |J with respect to the expected value of 0.0 is 
76.3 for 33 degrees of freedom. Removing sample number 
10 reduces the \ 2 to 57.0, and further removing sample 
number 11 reduces it to 48.7, which is just 2.25 standard 



deviations away from the expected value of 31.0. Therefore 
to the zeroth order of accuracy it is reasonable to assume 
that n 2 is the same (i.e., 0.0) for all samples in the in- 
tegrated profile. Its standard deviation is 0.0024 in the 
on-pulse window; then the standard deviation of \x can 
be estimated as ^0.0024/2 « 3.5% (by simple algebra). 
Thus the rms x-ray flux variation at any phase in the on- 
pulse window of the integrated profile of Crab pulsar is 
rs 3.5% of its mean value. Then the rms variation of the 
total on-pulse x-ray flux of the Crab pulsar will be ~ 3.5 / 
V^26 0.7%. A similar calculation for the entire integrated 
profile gives •\/0.0023/2/v / 33 w 0.6%, which is essentially 
the same result. This is a much tighter limit than the « 
7% quoted by Patt el al ( 199S ). These authors used to- 
tally 105 000 periods and a different method of analysis, 
on account of which their result might be dominated by 
photon noise. The result of this section should ideally re- 
flect the actual x-ray intensity variations of Crab pulsar 
(wave noise), that contain information about the physics 
of the x-ray emission mechanism. 

To the next order of accuracy appears to be cor- 
related with the shape of the integrated profile; both the 
lower panel of Fig. as well as the x 2 discussed above 
point to this. The fjr at the phase of the second peak of 
the integrated profile also appears to be enhanced. The 
on-pulse and off-pulse x 2 are 69.0 and 7.3, for 26 and 7 
degrees of freedom respectively; the former is 6 standard 
deviations away from the expected value, which is quite 
high. However more data, or better analysis, is needed to 
confirm this with good statistical significance. 

3.3. Giant pulses 

Comparison of photon counts in the off- and on-pulse win- 
dows of Fig. |l| shows that the Crab pulsar does not emit 
"giant pulses" at x-ray energies. The mean on-pulse win- 
dow photon rate is 17.57 photons in 26 time samples, 
while the maximum is about 31 photons. This implies 
that at x-ray energies the Crab pulsar emission occasion- 
ally increases by « 31.0 / 17.57 ~ 1.76 times, at the 
most, whereas at radio wavelengths the giant pulse en- 
ergy is about 10 to 100 times its mean value (Lundgren et 
al 1995). Further, the on-pulse photon distribution fits a 



Poisson distribution very well, and there is no discernible 
excess probability at higher pho ton r ates. Finally, fol- 
lowing the method of Ritchings ( 1976 ) and Vivekanand 
( 1995 ), a deconvolved photon distribution was obtained 
that represents the true on-pulse photon distribution of 
the Crab pulsar; it also does not indicate the presence of 
an excess probability at higher photon rates. 



4. Correlation of flux in the two peaks 



Patt et al (1999) report that the x-ray flux in the two 
peaks of the integrated profile of Fig. [y are uncorrelated. 
Their integrated profile consists of 16 samples, in which 
the two peaks are defined as covering phases 1 to 4 plus 
15 and 16, and phases 5 to 11, respectively, (Fig. 1 of 
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Patt et al 1999). They obtain the covariance between the 
number of photons in these two peaks to be —7.7 x 10 -4 , 
which they claim is consistent with zero. In this work the 
dotted curves of Fig. [l] are the models for the two peaks. 
The x-ray flux in each period was multiplied by the two 
curves, to obtain the flux in the two peaks, respectively. 
A normalized correlation coefficient (in contrast to the 
covariance of Patt et al 1999Q was computed for these two 
fluxes, across all 1868112 periods; the result is 3.5 (± 
10.0) x 10 -4 , which is consistent with zero. Therefore it 
is concluded that the x-ray fluxes of Crab pulsar in the two 
peaks of its integrated profile are indeed uncorrelated. 

5. Discussion 

The x-ray flux variations of Crab pulsar have a standard 
deviation of « 0.6% to 0.8% of its mean value over time 
scales ranging from a period to almost a month; that this 
number is similar over such wide time scales may or may 
not be a coincidence. This is consistent with the Crab pul- 
sar's behavior at optical and UV energies (Percival et al 



1993) and at IR energies (Lundgren et al 1995). Therefore 



the Crab pulsar flux variations at higher energies are in- 
significant compared to those at radio wavelengths, where 
the modulation index /j, is « 1.0, most of which apparently 
comes from the giant pulses (Lundgren et al 1995). Now, 



this causes problems for the conjecture that both the radio 
and the high energy emission are related through a com- 
mon electric current. If the flux variations of Crab pulsar 
are due to temporal variation in the number of basic emit- 
ters, then the radio flux variations would also have been 
at the ~ 0.7% level. If they are due to temporal variation 
in the coherence of the basic emitters, then the radio flux 
variations would have been at the 2 x 0.7 ~ 1.4% level, 
since the intensity of a coherent emission mechanism is 
proportional to the square of the number of basic emit- 
ters. If they are due to temporal variati on in the angle of 
some elementary beams (Lundgren et al 1995 ), one would 
have expected the radio flux variations to be smaller than 
those at high energy, since the common beams would be 
larger at the larger wavelengths. One could surely postu- 
late variations in the angle of radio beams only, but that 
would have to be justified on the basis of some other inde- 
pendent physical mechanism. In summary the difference 
in the radio and high energy flux variations of Crab pulsar 
is difficult to explain, if the basic charged emitters at both 
wavelengths are somehow related. 

Another method, of amplifying the very small flux vari- 
ations of Crab pulsar at high energies to the very large 
variations at radio wavelengths, would be to somehow 
use a fraction of the « 10 7 amplification factor of par- 
ticles in the gaps, due to cascading e + -e _ pair produc- 



tion (Ruderman & Sutherland 1975; Cheng et al 1986a 



1986b). One could probably postulate that the radiation 



that is emitted by the later generation of charges in the 
pair cascade process suffers greater variation in its inten- 
sity, due to amplification of the variation in the number 
density of charges. This might also imply that one should 



see a monotonic increase in the modulation index as one 
observes at larger wavelengths. More detailed study of 
Cheng et al (1986a, 1986t) and Romani & Yadigaroglu 
(1995) models is required for a reasonable solution. In any 
case, the results of this paper set very strong constraints 
on the explanation for the relative flux variations at the 
radio and x-ray energies. 

The possible correlation of /i with the integrated pro- 
file in Fig. |4|, if confirmed in future, might set very strong 
constraints on the basic emission mechanism of high en- 
ergy emission from the Cra b pulsar. In the framework of 
Cheng et al ( 1986a , 19861 ) and Romani & Yadigaroglu 
( 1995| ) models, the two peaks of the integrated profile are 
cusps created by emission from different magnetic field 
lines, that add in phase along different directions due to 
relativistic aberration. A well defined relation between the 
mean flux and its variance at each point in the integrated 
profile (for example, /i might vary as (/)") will be an ad- 
ditional constraint, along with the exact shape (I) of the 
integrated profile, on the x-ray emission mechanism. 

Further work on these data is in progress, that studies 
issues such as (a) verifying if the Crab pulsar shows at 
x-ray energies the three phenomenon that are often seen 
in several radio pulsars - "pulse nulling" "systematic sub 
pulse drifting" and "mode changing" ; (b) looking for spe- 
cial behaviour in the x-ray integrated profile at the phase 
of the radio precursor, which is supposed to be different 
from the rest of the radio integrated profile; (c) compari- 
son of the peak and the bridge x-ray emission of the Crab 
pulsar, which might further constrain the models of Cheng 
et al ( |T986a| [1986% Romani & Yadigaroglu fll995| ), and 
Cheng et al p000| . 



On the theoretical front, it is probably worth exploring 
the simultaneous modeling of (/) and aj (or /i) at x-rays 
for rotation powered pulsars, and more specifically for the 
Crab pulsar. 
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